Double transverse spin asymmetry in the p^p^ Drell-Yan process from Sivers functions 
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We show that the transverse double spin asymmetry (DSA) in the Drell-Yan process contributed 
only from the Sivers functions can be picked out by the weighting function jfe (cos(0 — cj>s 1 ) cos(<^> — 



6g 2 ) + 3sin( 



) sin(<ji> — <t>s 2 ))- The asymmetry is proportional to the product of two Sivers 



functions from each hadron /£y x fy^. Using two sets of Sivers functions extracted from the 
semi- inclusive deeply elastic scattering data at HERMES, we estimate this asymmetry in the 
Drell-Yan process which is possible to be performed in HESR at GSI. The prediction of DSA in the 
Drell-Yan process contributed by the function gir(x, k§d, which can be extracted by the weighting 
function ^p-(3cos(<^ — 4>Si) cos(0 — 4>s 2 ) + sin(0 — <ps t ) sin(^) — <^s 2 ))i is also given at GSI. 

PACS numbers: 13.88.+e, 13.85.Qk 



I. INTRODUCTION 

The Sivers effect [l[ was proposed originally to explain 
the large single spin asymmetries (SSA) observed in in- 
clusive pion hadro-production (p^p — > irX) at FNAL 
The effect can be quantitatively described by a k-r- 
dependent distribution named as Sivers function [E H[ 
/j^(x,k|,), which is the distribution of unpolarized par- 
tons in a transversely polarized proton. It arises from 
a non-trivial correlation between the nucleon transverse 
spin and the intrinsic transverse momenta in the nucleon. 
Despite its (naively) T-odd property H , Sivers function 
has been proven to be non- vanishing [6| due to its special 
gauge -link property 0,1,11. 

Recently the SSA measured in semi-inclusive deeply 
inelastic scattering (SIDIS) processes with transversely 
polarized targets at HERMES 0,0 and COM- 
PASS [H, [lj], has been shown to be interpreted by 
the Sivers effect. The asymmetry is identified by the 
angular dependence sin(</> — cj>s), where and <ps de- 
note respectively the azimuthal angles of the produced 
hadron and of the nucleon spin polarization, with respect 
to the lepton scattering plane. The coexistent Collins 
asymmetry [f|, with a angular dependence sin(</> + 4>s), 
has also been measured in those experiments. The 
data on the Sivers SSA has been utilized by different 
groups 
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to extract the Sivers functions 
of the proton, especially those for the u and d quarks, on 
the basis of the generalized factorization [2(| Hl[ . Those 
sets of parametrization of the Sivers functions are qual- 
itatively in agreement [22| among themselves, and were 
applied to predict the Sivers SSA in various processes in 
the established or planed facilities, such as the SIDIS at 
JLab, and the Drell-Yan process at COMPASS, RHIC 



and GSI. 

In this paper, we will investigate the role of the Sivers 
function on the transverse double spin asymmetry (DSA) 
in the Drell-Yan process. The transverse DSA has been 
investigated [23| for many years, and is believed to have 
advantage to unravel the transverse spin property of 
the nucleon [24J, especially the transversity distribution 
hi(x) [2o| . Various azimuthal asymmetries contributed 
by different kf-dependent distribution functions have 
been analyzed and given in Refs. [IE] and [271 ] . As 
shown in Ref. 27], the Sivers function contributes to 
the DSA in the Drell-Yan process through the product 
fir x /it- However, this DSA is mixed with the contri- 
bution from another k^-dependent distribution function 
gix(x,]i^). We will show that through the appropri- 

ate weighting function ^-(cos(0 — cos(</> — </>s 2 ) + 
3sin(0 — 4>Si) sin(0 — 4>s 2 ))i the- asymmetry from the 
Sivers function can be isolated without mixing with the 
contribution from other functions. Using two sets of pa- 
rameterizations [IE HH of the Sivers functions we cal- 
culate the double spin asymmetry from the Sivers func- 
tions in the p^p^ Drell-Yan process at GSI. An asymme- 
try around 1 % is predicted. The asymmetries estimated 
from these two sets of Sivers functions are quantitatively 
different. Therefore measuring the DSA in the Drell-Yan 
process can provide new information on Sivers functions, 
especially their sizes. The transverse DSA contributed 
by fllT^! kfi) through the product giT x grr can also 
be picked out by another weighting function. We esti- 
mate this asymmetry by adopting a gvr coming from the 
combination of a Lorentz invariance relation presented 
in Refs. [H, [29j and the Wandura-Wilzeck approxima- 
tion [ai. 



II. EXTRACTING DSA CONTRIBUTED BY 
THE SIVERS FUNCTIONS 
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The importance of the transverse-momentum distri- 
butions of quarks for a full understanding of the struc- 
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ture of hadrons has been widely recognized in the last 
decade ^H^HEUll- A comprehensive leading- twist tree 
level analysis of the (spin dependent) Drell-Yan process 
in terms of k^-dependent distributions has been given 
in Ref. [26[ . The role of the T-odd k^-dependent distri- 
butions in this process has been presented in Ref. [27l |. 
In the Collins-Soper frame [33| the leading order unpo- 
larized differential cross-section for the Drell-Yan process 

h i( p i) + MP2) -> 7*(?) + x -» + + X has 

the form [2711 
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+B(y)cos2(f)T 



(2h ■ prh • ky — Pt ■ k^ 



hf q hf 9 
MiM 2 



(1) 



where q denotes the quark flavors, the notation 

T[hh] = J dV^k^^pT + kT-qr) 

x/iOcx.PiOAte.kr) (2) 

shows the convolution of transverse momenta, Q 2 = q 2 
is the invariance mass of the lepton pair, qx is the trans- 
verse momentum of the lepton pair, h = / Qt, <j> is the 
angle between the hadron plane and the lepton plane, 
and 



•U/;i = ( \ -y + y 2 ) = J(i + cos 2 i 



B(y) = y(l-y) 
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(3) 
(4) 



in the cm. frame of the lepton pair. 

The function hj_in the second line of |T]) is the Boer- 
Mulders function [J] , the chiral-odd partner of the Sivers 
function. This function has attracted a lot of interest [34|, 
IH, [36| recently because it can account for the anomalous 
cos 2cj) asymmetries [37| H|| observed in the unpolarized 
Drell-Yan process, given by the second term of Eq. ([1]) 
has shown. 

The leading order differential cross-section for the dou- 
ble transversely polarized Drell-Yan process is [l?) 



da {2 \h\h\ -> UX) a 2 
dfldxidx2d 2 q T 

, Mv) 



3Q 



2 



E 



|S 1T | |S 2T | cos(20-0 Sl -<As 2 ) 



fXq fJ-q _ q q 
J IT J IT yiTi/lT 



h-pT h k T 



MiM 2 



Ax{y) 



SitI Pari 



cos 



<AsrJ cos(0- 0s 2 ) 



f ±q r-Lq 
JIT JIT 

MiM 2 



Mv) 



'IT I 



3 2T 



sm 



0sJ sin(<^- 0s 2 ) 



9 g 

9\t9\t 

MiM 2 



(5) 



The . . . indicates the terms which will not contribute in 
our analysis below, <pg 1 and (f>s 2 are the angles between 
^it 1 S2T an d the lepton plane, respectively. 

As shown in ([5]), the Sivers function can contribute 
to the transverse DSA through the product x f^ T . 
However this asymmetry is mixed with the asymmetry to 
which it contributes another k^-dependent distribution 
gyrix, kfi). The main goal of this paper is to isolate 
the asymmetry contributed by the Sivers function. The 
starting point is the method introduced in Ref. [3^ |. by 
which one integrates the differential cross section with a 
proper weighting function W(Qt, 4>, <f>Si , 4>S 2 )i as follows: 

(W{QtA^ Si As 2 )) 



4 1 4>Si qr (ifldxidx2d 2 q T 
xW(Qt,<£,0Si,0S 2 )- 



(6) 



With the above weighting procedure, one can pick up the 
terms in which one is interested. Besides this, one can 
de-convolute the transverse momentum integration in a 
model independent way. 

The unpolarized angular independent cross section can 
be picked out by using the weighting function 1, from 
Eq. Q: 
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We denote Wc 
$m{4>- <j) Sl )sm.(<) 



cos (9 



- 4> Sl ) cos(0 - <j) S2 ) and W s = 
Given the weighting function 

u , .. ! ,...•••....».■•,, .... 1 *.,l2 — M, i.e. the colliding two 
hadrons are nucleons), we can obtain the following term 
from (O): 



o 
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where fix (x) and g\^ (x) are the first k^-moments, de- 
fined as: 

At (1) (*) = fd^T^fM^r), (9) 
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g$(x) = d'k 



' 2M 2 



gi T (x,\s%.). 



(10) 



The factor Q T introduced in the weighting function en- 
sures that the transverse momentum integration in ([8]) 
can be de-convoluted (for details, refer to the Appendix). 

Again, applying the weighting function j&Ws on ([5]), we 
arrive at 

( A(yW em y\/Ql \ 

= -^E e '{ 3 [At (1)9 (^i)/it (1)? (^)- 
q 

^\ Xl ) 9 {^(x 2 )]-2g{^( Xl )g[^(x 2 )} 



TT 



At 



g^ 9 ( X i)g\^( X2 ) 
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(11) 



Therefore, combining ([8]) and (jlip. we can extract the 
term contributing to the transverse DSA and coming only 
from the Sivers functions: 



V 3Q 2 



TT 



-8^^e 2 q f^"(x 1 )f^(x 2 ), 



(12) 



with the weighting function jp(Wc + 3Ws). 

By taking the ratio between (fl~3|) and 0, we define 
the weighted double spin asymmetry as follows 
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(13) 



E,^/i 9 (xi)/i*(x 2 ) 

The above equation thus provides a possibility to study 
the Sivers function by measuring the transverse DSA in 
the Droll- Yan process. 

Also, from ©, (fTTj) and (J7J) we can get another type 
of DSA: 



^(iWc + Ws) 



TT 



Win 



which is contributed only by gi^. 



(14) 



III. NUMERICAL RESULTS 

In this section we will give numerical results on the 
DSA from the Sivers functions. We consider the trans- 
versely polarized proton antiproton Drell-Yan process, 
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FIG. 1: The DSA in the proton antiproton Drell-Yan process 
at GSI coming only from the Sivers functions, and calculated 
from Eq. ([13]). The kinematics is s = 45 GeV 2 and Q 2 = 
2.5 GeV 2 . The solid and dashed curve use the Sivers functions 
in Ref. [l6l ] and in Ref. [jq|. respectively. 



where the valence Sivers functions are involved, so that 
a larger asymmetry should be measured compared to 
the p^pT Drell-Yan process. The p^p^ Drell-Yan process 
is possible to be performed in the planned high energy 
storage ring (HESR) [U at GSI. We study the trans- 
verse DSA at GSI from the Sivers functions, based on 
Eq. (fTBf . For this end we need to know the input for 
the Sivers functions. Several groups [l(| [H, [l^] have 
parameterized the Sivers functions based on the data of 
SIDIS at HERMES [HG3, and partially based on COM- 
PASS data [13j . The kinematics in GSI can be chosen as 
the cm. energy s = 45 GeV 2 . For the invariance mass 
square of the lepton pair we choose Q 2 = 2.5 GeV 2 , which 
is close to the scale at HERMES. Therefore these sets of 
Sivers functions extracted from the data of HERMES can 
be applied to predict the asymmetries at GSI in the kine- 
matics regime we give above. We will adopt two sets of 
Sivers functions, which are the sets in Ref s. fl6ll and [l||, 
respectively. The Sivers functions in Ref. [19| can not be 

applied here since in that paper / r ^ 1 ^ 2 ' ) [x) is given while 
we use f^f{x) in our calculation. 

To use these Sivers functions one should notice that 
T-odd distribution functions in the DIS and in the Drell- 
Yan process have a minus sign difference 0] ■ However in 
the p^p^ Drell-Yan process two Sivers functions appear in 
the product, therefore the sign difference doesn't matter 
here and the functions can be used directly. 

In Ref. the Sivers functions are parameterized as 



f 1T q ( X , k 2 -) 



N q ( x )r i (x)gCk T )H^T), 
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extracted compared to the parametrization in Eq. (|19p . 
which will lead a larger asymmetry compared to the 
dashed curve in Fig. [TJ Thus the difference between the 
asymmetries from the two sets of Sivers functions may be 
reduced. Depending on the accuracy of the experimental 
measurements on the transverse DSA at GSI, useful 
constraints on the Sivers functions could be obtained, 
but it might be hard to distinguish between different 
parameterizations without high precision measurements. 

Finally, we will predict the DSA contributed by the 
function g\x{x, k 2 ^) at GSI. This function, describing lon- 
gitudinal polarization of quarks in the transversely po- 
larized target, also plays role in the double polarized 
(longitudinal-transverse) SIDIS process [H, |43[ . A treat- 
ment on <7iT(a;,k^) is the so-called Lorentz invariance 
relation that connect the first k^ moment of giT{x,\i^) 
with the twist-three distribution function g% (x): 



FIG. 2: The DSA in the proton antiproton Drell-Yan process 
at GSI coming from the function gir, defined in Eq. (I14|l . 
The kinematics is s = 45 GeV 2 and Q 2 = 2.5 GeV 2 . 



with 



N q (x) = N q x a " (1 - x) b " 



(a q + b q y a *+ b ^ 



dq O q q 



(16) 
(17) 



for q = u,d. For the function /i(k^) two options are 
considered: 



(a) h(l4-) 



2k T M 
k% + Ml 



(b) VTe^e-^'. (18) 



In our calculation we will adopt option (b) in Eq. (fT5|) . 
and the central values of their fit. This parametriza- 
tion has taken advantage of the more precise data [l2j at 
HERMES. 

In Ref. [18] the authors give the set of Sivers functions 
for the u and d quark as, 

z/it (1) '» = -xfiT 1] ' d = -0.17x a66 (l - x) 5 , (19) 

extracted from the published HERMES data [III) and 
whose form is based on the limit of a large number of 
colours N c . 

For the unpolarized distribution we use the 
MRST2001(LO set) parametrization In Fig. Q] 

we present the DSA from Sivers functions at GSI, as a 
function of x\. A sizable asymmetry is predicted. The 
asymmetry (solid line) based on the Sivers functions from 
Ref. [IH is much larger than the asymmetry (dashed 
line) based on the Sivers functions from Ref. [l8| . As ex- 
plained in Ref. [42j , taking into account the more precise 
data [12| of HERMES, larger Sivers functions can be 



(20) 



Using the Wandzura and Wilczek approximation for g\ 
g q 2 {x)^-g a 1 {x)+ f dy?M, (21) 



y 



the following relation was derived in Ref. 
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(22) 



For the polarized parton distribution we apply the 
GRSV2001 (standard scenario) parametrization [4J|, and 
for the unpolarized distribution we use GRV98 LO 
parametrization [45j . following the choice in Ref. [43| . 
In Fig. [5] we show the DSA contributed by grr{x, ky) in 
the p1p1 Drell-Yan process at GSI with s = 45 GeV 2 and 
Q 2 = 2.5 GeV 2 . An asymmetry of 1% is predicted. 

We end this section with some comment. In our cal- 
culation, especially in the case of Siver DSA, we choose 
Q 2 = 2.5 GeV 2 . This value is consistent with the aver- 
aged scale (Q 2 ) in the HERMES experiment, from which 
the Sivers functions were extracted. Therefore, the pa- 
rameterizations for Sivers functions in Refs. [3, can 
be applied here without further assumptions. Experi- 
mental measurements at GSI can also cover the contin- 
uous Drell-Yan masses 2 — 5 GeV which corresponds to 
Q 2 in the range 4 — 25 GeV 2 . To estimate the asym- 
metries in this region one should use the fitted functions 
evolved to the relevant scale, which is not trivial for the 
k^-dependent distributions [46]. Therefore we assume 
that the ratios in Eqs. (fT"3|) and (fl"4]) scale with Q 2 . In 
this region, The result is similar to the one which can be 
obtained at the fixed value of Q 2 = 2.5 GeV 2 . Also there 
is the subtlety that the next to leading order correction of 
the hard process could lead the substantial JGfactor on 
the transversely polarized cross-section. Since we calcu- 
late an asymmetry, which is essentially a ratio where the 
Q 2 dependences in the numerator and denominator tend 
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to cancel each other, the effect of both the Q 2 dependence 
and K-factors do not introduce a strong influence on the 
resulting prediction coming from Eqs. (fT3|) and (fl~4|) . 



IV. SUMMARY 

We have performed an analysis of the transverse DSA 
in the Drell-Yan process contributed by the Sivers func- 
tions through the term f^p x f^ T . The asymmetry can 
be isolated through the appropriate weighting function 

ffc (cos(<j)-(j) Sl ) cos(0- 4>s 2 ) +3 sin(0-(/) Sl ) sm(<£- </>s 2 )), 
without mixing with the contribution from other distri- 
bution functions. Using two sets of Sivers functions pa- 
rameterizing the SSA data in the SIDIS process, we cal- 
culate the double spin asymmetry in the p^p^ Drell-Yan 
process from the Sivers functions at GSI. An asymmetry 
around to 1 % is predicted. The asymmetries estimated 
from these two sets of Sivers functions are quantitatively 
different. Therefore measurements of the DSA in Drell- 
Yan process can provide new information on the Sivers 
functions, especially their sizes. The transverse DSA con- 
tributed by giT(x,k T ) through the product gyr x 9it in 
the Drell-Yan process can also be picked out by a weight- 
ing function. We estimate this asymmetry at GSI by 
adopting gyr from the combination of the Lorentz invari- 
ance relation and the Wandura-Wilzeck approximation. 
The investigation on the double transversely polarized 
Drell-Yan process thus can shed light on the knowledge of 
ky-dependent distribution functions, including the Sivers 
functions. 
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APPENDIX: MOMENTS 

To derive ([5]) and (fTTj) we have used the following trans- 
verse momentum integrations: 

f Q 2 

I d 2 k T d 2 p T <5 2 (q T - k T - p T )j^(k T • p T ) 

xf{x 1 ,k T )f(x 2 ,p 2 T ) 

J d 2 k T d 2 p T (k T + p T ) 2 k T • Pt 



M 2 

<} 
2 



x/(xi,k|)/(x 2 ,p|) 

i/2 Jd 2 k T d 2 p T (k T ■p T ) 2 f(x 1 ,k 2 T )f{x 2 ,p 2 T ) 
~M 2 _ 

x/(a; 1 ,k|)/(x2,Px) 



//•')•') 2 2 

d 2 k T d 2 p T [kip p T + k T p T 

C!,k|)/(x2,Px) 
4M^«(x 1 )/( 1 )(x 2 ). (A.l) 



f Q 2 - 

/ d 2 k T d 2 p T 5 2 ((\ T -k T ~ Pt) • k T h • p T 

xf(x u k 2 T )f{x2,V 2 T ) 
= J d 2 k T d 2 p T (k T + p T ) ■ k T (k T + p T ) ■ p T 

xf(x 1 ,k 2 T )f(x2,p 2 T ) 

= W I d2k T d2 Pr(p T K + (kr • Pt) 2 ) 

xf(x 1 ,k T )f(x2,PT) 
= 6M 2 / (1) (*i)/ (1) (*2) (A.2) 

In the above integrals, the terms containing odd numbers 
of k T or p T vanish after being integrated over k^ or p^- 
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